Spin filling of valley-orbit states in a silicon quantum dot 
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We report the demonstration of a low-disorder silicon metal-oxide-semiconductor (Si MOS) quan- 
tum dot containing a tunable number of electrons from zero to N = 27. The observed evolution 
of addition energies with parallel magnetic field reveals the spin filling of electrons into valley-orbit 
states. We find a splitting of 0.10 meV between the ground and first excited states, consistent with 
theory and placing a lower bound on the valley splitting. Our results provide optimism for the 
realization in the near future of spin qubits based on silicon quantum dots. 
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1. Introduction 

Semiconductor quantum dots pQ are islands to which 
electrons can be added one by one by means of an elec- 
tric field. Like real atoms they have discrete quantum 
levels and can exhibit phenomena such as shell filling [2] , 
where orbital levels are filled by spin-paired electrons to 
produce a spin-zero many-electron state. Quantum dots 
also provide a promising platform for spin qubits, which 
can have long coherence times due to the weak coupling of 
spins to local fluctuations in charge. For a quantum dot 
to be useful as a spin qubit it is essential to understand 
the details of its excitation spectrum and its spin-filling 
structure. One powerful method to probe the spin fill- 
ing is via magnetospectroscopy. This has been applied 
to both vertical [2] and lateral GaAs/AlGaAs quantum 
dots [3 , showing ground-state spin filling in agreement 
with Hund's rule. 

Silicon devices are attractive for spin-based quantum 
computing [HE] and spintronics [6] because of their scala- 
bility and long spin coherence times [7 . Silicon quantum 
dots, in particular, have potential as electron-spin qubits, 
but to date it has not been possible to create devices with 
the low disorder present in their GaAs counterparts [2j [3] . 
This is primarily due to disorder at the Si/Si02 interface, 
which has made it difficult to achieve single-electron oc- 
cupancy [HIS]. 

In addition, the conduction band structure in silicon 
is complex and only a few experiments have been car- 
ried out to examine the spin states in either Si MOS or 
Si/SiGe quantum dots [SHT3], The valley degree of free- 
dom makes the measurement and interpretation of spin 
states in all silicon-based dots non-trivial p~4j [15] , while 
for Si MOS dots the substantial amount of disorder usu- 
ally present at the Si/Si02 interface impedes the ability 
to make smooth potential wells. 

In this work we present the investigation of a Si MOS 
quantum dot with lower disorder than any studied to 
date, in which it is possible to analyse the electron oc- 
cupancy in a manner previously inaccessible. We deduce 
the spin filling of the first 12 electrons in the dot from 



ground- state magnetospectroscopy measurements. The 
formation of a two-electron (N = 2) spin-singlet state at 
low magnetic fields confirms that there is no valley de- 
generacy present, while the magnetic field dependence of 
the higher-order Coulomb peaks allows us to deduce the 
level structure for the first four electrons. 

In the following section we present the architecture 
of the quantum dot and the charge stability diagram in 
the few-electron regime. We then, on section 3, study 
the valley-orbit states in this quantum dot and extract 
a valley-orbit splitting of 0.10 meV. In section 4 we 
investigate the spin filling of the first 4 electrons in this 
quantum dot in detail. We then present the spin filling 
of the 5 th to 12 th electrons, discussing some anomalies 
observed, before concluding in section 6. 

2. Low-disorder silicon MOS quantum dot 

The triple-layer gate stack in our structure (Figure la 
and lb) provides excellent flexibility for tuning the bar- 
rier transparency and the energy levels of the dot in- 
dependently, see Supplementary Information for fabrica- 
tion processes. The lowest layer defines the barrier gates 
(Bl and B2). They are used to define the dot spatially 
and control the tunnel coupling. The second layer of 
gates defines the soure-drain leads (LI and L2). The 
lead gates induce the electron accumulation layers that 
act as source-drain reservoirs. The plunger gate (P) ex- 
tends over the barrier gates, lead gates and the dot is- 
land, and is used to control the electron occupancy of the 
dot. Figure lc is a plot of the differential conductance 
dl/dVsB of the device versus plunger gate voltage Vp and 
source-drain voltage Vsd, showing the familiar "Coulomb 
diamond" charge stability map. Before the first charge 
transition the diamond edges open entirely to a source- 
drain voltage |Vsd| > 20 mV, because the quantum dot 
has been fully depleted of electrons. We have previously 
reported a device with similar gate architecture but an 
accidental parallel quantum dot created distortion of the 
charge stability map in the few-electron regime, compli- 
cating the interpretation of the dot's level structure [8]. 
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FIG. 1: a, Scanning electron microscope image and b, 
Schematic cross-section of a Si MOS quantum dot. c, Sta- 
bility diagram of the device in the few-electron regime. By 
decreasing the plunger gate voltage Vp , electrons are depleted 
one-by-one from the dot. The first diamond opens up com- 
pletely indicating that the first electron has tunneled off the 
dot. d, Coulomb oscillations as a function of plunger gate 
voltage Vp for the first 27 electrons in the dot. Vp is com- 
pensated by Vb2 to suppress the non-monotonic background 
conductance. Inset: Addition voltage V a dd versus electron 
number N calculated as the difference between two consecu- 
tive Coulomb peaks in plunger gate voltage. 



Here, clear and sharp Coulomb peaks mark the first 27 
electrons entering the dot, see Fig. Id, while the charge 
stability map of Fig. lc shows no distortions from disor- 
der potentials. 

As with quantum dots in GaAs/AlGaAs [2], shell 
filling has very recently been observed in Si/SiGe 
quantum dots, with a filled shell structure observed for 
N = 4 electrons [12 . The addition spectrum of our 
Si MOS quantum dot (inset of Figure Id) also shows 
a noticeable peak at N = 4. A filled shell at TV = 4 
would be consistent with the filling of a first orbital 
state in a two-valley system, however, an accurate 
description of orbital and valley levels in silicon quan- 
tum dots is somewhat more complex, as described below. 

3. Valley-orbit splitting 

In recent years, valley physics in silicon has been stu- 
died extensively both theoretically [15-22 and experi- 
mentally [23 - 26 . In bulk silicon, there are six degenerate 
conduction band minima (valleys) in the Brillouin zone, 
as depicted in Figure 2a. Confinement of electrons in 
the z-direction at the Si/Si02 interface lifts the six- fold 
valley degeneracy: four A- valleys with a heavy effective 
mass parallel to the interface have an energy several tens 
of me V higher than the two T- valleys [27] . The sharp and 



fiat interface produces a potential step in the z-direction 
and lifts the degeneracy of the T-valleys in two levels 
separated by the valley splitting Ey. Theoretical predic- 
tions for the valley splitting are generally on the order 
of 0.1-0.3 meV [15 , 20 j. Experimental values in Si inver- 
sion layers mostly vary from 0.3-1.2 meV [23 . A very 
large valley splitting of 23 meV in a similar structure has 
also been measured [24] and is explained in [20] [21] . Re- 
cently, resonant tunneling features spaced by ~0.1 meV 
in a single-crystal silicon quantum dot were attributed to 
valley excited states [26], while measurements on Si/SiGe 
quantum dots revealed valley splittings in the range of 
0.12-0.27 meV [12] . 

Valleys and orbits can also hybridise [19], making it 
inappropriate to define distinct orbital and valley quan- 
tum numbers. Depending on the degree of mixing, the 
valley-orbit levels behave mostly like valleys or like or- 
bits. Instead of referring to a pure valley splitting we 
therefore adopt the term valley-orbit splitting, AEyo = 
Eyo2 — Eyoi for the difference in energy between the 
first two single-particle levels, Even and Evc>2- This is 
sometimes referred to as the ground-state gap [19 . 

Full electrostatic control of the electron number allows 
us to investigate the spin filling by measuring the mag- 
netic field dependence of the electrochemical potential 
/iTV, which is by definition the energy required for adding 
the N th electron to the dot. The slope of /jLjsf(B) is given 
by [SI 

^ = -gn B AS tot (N), (1) 

where g is the g-factor, the Bohr magneton /iB = 58 
fieV/T and AS tot (N) is the change in total spin of the 
dot when the N th electron is added. The electrochemical 
potential has a slope of -fg/ie/2 when a spin-up electron 
is added, whereas addition of a spin-down electron results 
in a slope of —g^/2. The rate at which [in changes with 
magnetic field thus reveals the sign of the added spin. For 
the experiments in this work we apply the magnetic field 
B parallel to the Si/Si02 interface. 

The conductance at the first two charge transitions 
is plotted as a function of the electrochemical potential 
energy and the magnetic field in Figure 2b. Here, the 
Coulomb peak positions in gate voltage are converted to 
electrochemical potential [In using the lever arm ap ex- 
tracted from the corresponding Coulomb diamonds. The 
blue lines above the Coulomb peaks are guides for the 
eye with slopes of ig/ie/2, as predicted by equation (1) 
using g = 2 for bulk silicon. Since the first Coulomb 
peak moves down in energy with increasing magnetic 
field the peak corresponds to a spin-down electron enter- 
ing the quantum dot, as expected for the N = 1 ground 
state. For B > 1 T the second Coulomb peak also falls 
in energy with increasing B at a rate close to — g/iB/2, 
however, for low magnetic fields the peak noticeably in- 
creases in energy with £?, leading to a "kink" (marked 
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FIG. 2: a, Conduction band minima (valleys) in bulk silicon, 
showing six ellipsoids of equal energy in the Brillouin zone. 
Each ellipsoid has two light traverse mass (mt) and a heavy 
longitudinal mass (mi). Under the z-direction confinement at 
the Si/Si02 interface, the six- fold degenerate valleys split into 
two T- valleys (lower in energy) and four A-valleys (higher in 
energy). The sharp interface potentials split the T- valleys by 
an amount Ey . b, Magnetospectroscopy of the first two elec- 
trons entering the quantum dot. The circle 2a marks a kink in 
the second Coulomb peak at ^0.86 T. The arrows in the boxes 
(VOi for valley-orbit 1 and VO2 for valley-orbit 2) represent 
the spin filling of electrons in the quantum dot. Coulomb peak 
positions in gate voltage are converted to energies using the 
lever arm ap extracted from the corresponding Coulomb dia- 
monds, c, A model showing that the valley-orbit splitting can 
be estimated from the magnetic field at which AEyo = AEz , 
i.e. when the spin-up state of VOi is at the same energy as 
the spin-down state of VO2. For B < 0.86 T, the first two 
electrons fill with opposite spins in the same valley-orbit level 
(left panel). As we increase the magnetic field, the Zeeman 
energy exceeds the valley-orbit splitting and the second elec- 
tron occupies a spin-down state in valley-orbit 2. The sign 
change appears as a kink and occurs when the valley-orbit 
splitting is equal to the Zeeman energy (0.10 meV). 



netic field is applied, the spin-down and spin-up states 
are split by the Zeeman energy Ez- Above 0.86 T the 
spin- up state of valley-orbit level 1 (VOi) is higher in 
energy than the spin-down state of valley-orbit level 2 
(VO2) and it becomes energetically favoured for the sec- 
ond electron to occupy the latter, i.e. VO2. At the 
kink the valley-orbit splitting equals the Zeeman energy, 
which is 0.10 meV at 0.86 T. With the interfacial elec- 
tric field of ~2xlO T V/m extracted from Technology 
Computer-Aided-Design modeling for our device struc- 
ture, the valley-orbit splitting agrees well with modeling 
results (0.08-0.11 meV) based on the effective- mass ap- 
proximation [20 , 22 . We note that if no valley-orbit mix- 
ing were present, then A£Vo =0.10 meV would place a 
lower bound on the valley splitting for this structure. 

For B > 0.86 T the first two electrons fill two differ- 
ent levels split by AEyo = 0.10 meV. We note that the 
presence of a doubly degenerate ground-state level would 
demand the two electrons to exhibit parallel spin filling 
starting from T, since the two electrons would then 
occupy two different valley states in order to minimise 
the exchange energy [14 j. A valley- degenerate state is 
therefore ruled out by the results in Figure 2b. 

To assess the degree of valley-orbit mixing we compare 
the expected values for the orbital level spacing and the 
valley splitting. As stated above, theoretical calculations 
of the latter predict 0.1-0.3 meV. An estimate of the 
orbital level spacing in a quantum dot is given by 
27rh 2 /g v g s m*A pQ, where g v (g s ) is the valley (spin) 
degeneracy, m* the electron effective mass and A the 
dot area. For non-degenerate valleys, g v = 1 and g s = 2. 
Using the effective mass of 0.19mo, and the lithographic 
dot area of ~ 30 x 60 nm 2 we obtain an expected orbital 
level spacing of 0.7 meV. This value is considerably 
larger than the lower bound on the valley splitting, 
suggesting that the first two levels may be valley-like, 
however, to maintain generality we will continue to refer 
to the levels as valley-orbit states. 



2a) at B ~ 0.86 T. This kink (2a) is confirmed by sev- 
eral repeated measurements over positive and negative 
magnetic field, see Supplementary Information (Fig. S2). 
These results imply that at low magnetic field (before the 
kink), the second electron fills the quantum dot with its 
spin up. As we increase the magnetic field (after the 
kink), the sign of the second electron spin changes from 
up to down at B ~ 0.86 T. We note that in previous 
measurements on a similar quantum dot device, disorder 
and instability made it difficult to accurately probe this 
kink feature [8]. 

We explain the sign change observed here with a sim- 
ple model where the two lowest valley-orbit levels are 
separated by the valley-orbit splitting AEyo, see Figure 
2c. At zero magnetic field, the first two electrons fill 
with opposite spins in valley-orbit level 1. When a mag- 



4. Spin filling of the first 4 electrons 

We now turn to the spin filling for N > 2 electrons. 
Figure 3a shows the differential conductance as a func- 
tion of plunger gate voltage and barrier gate voltage Vb2 • 
The highly regular pattern of parallel Coulomb peak lines 
again demonstrates the low disorder in this device. In 
order to determine the spin filling for higher electron 
numbers we investigate the difference between successive 
electrochemical potentials as a function of magnetic field. 
The resulting addition energies E add (N) = /i/v — Hn-i 
have slopes which depend on the spin filling of two con- 
secutive electrons, according to [28] 
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FIG. 3: a, Differential conductance dI/dV s d as a function 
of barrier gate voltage Vb2 and plunger gate voltage Vp at 
B = 0T. The regular parallel Coulomb peaks are a signature 
of low disorder, b, Addition energies of the 3 rd and 4 th elec- 
trons versus magnetic field. Kinks are reproducible and ap- 
proximately symmetric over positive and negative magnetic 
fields, c, A simple model showing the evolution of single- 
particle energy levels Eyoi of valley-orbit i assuming only the 
Zeeman shift. Each level splits into two levels Evoi ± \g\i^>B 
at non-zero magnetic fields. The level crossings fit the kinks 
observed in the first four Coulomb peaks shown in Fig. 2b 
and Fig. 3b. 



where the first (second) arrow depicts the spin of the 
(N — l) th (7V th ) electron respectively. 

Figure 3b plots the measured addition energies, 
E&dd(N) = [In — Miv-i> for N = 2 to TV = 4 electrons for 
magnetic fields B in the range of —8 T< B < 8 T. We see 
that the data in Fig. 3b tend to follow dE a( ±a(N)/ 'dB = 
0,=b<7/iB, as expected from Equation (2). Furthermore, 
the E^a^N) data is relatively symmetric about = 0, 
indicating that the trends are real and not measurement 
artefacts. As a guide to the eye, we also show lines with 
slopes of exactly 0, ±g\iB (blue lines in Fig. 3b) that we 
interpret the E a ^(N) to be following. While in regions 
the match is not exact, we propose that these trend lines 
are the best qualitative fit to the data. We are thus able 
to infer spin states for each of the first 4 electrons at all 
values of magnetic field \B\ < 8 T. These spin states are 
labelled with red (green) arrows, representing spin down 



(up), in Fig. 3b. 

We now focus on the spin states of the these four elec- 
trons, N = 1 to N = 4. At low magnetic fields (< 0.8 T), 
the electrons populate the quantum dot ground states 
with alternating spin directions: |, t, |, t- Conversely, at 
high magnetic fields (> 4 T) a configuration with four 
spin-down electrons has least energy: |, |, |, |. Recently, 
parallel spin filling in a Si quantum dot was explained as 
a result of a large exchange energy and an unusually large 
valley splitting of 0.77 meV [9]. When the level spacing 
is smaller than the exchange energy, it is energetically 
favoured for two electrons to occupy two consecutive lev- 
els with the same spin sign. This is not the case for the 
device measured here: the anti-parallel spin filling of the 
first two electrons below 0.86 T is only possible in case 
of a small exchange energy (less than AEyo). This is 
an unexpected result for a dot of this size where the ex- 
change energy is predicted to be larger than the orbital 
level spacing [14]. Possibly the Coulomb interaction in 
the dot is strongly screened by the plunger gate. This is 
not unlikely since the distance from gate to dot (10 nm) 
is smaller than the dimensions of the dot itself (30—60 
nm). 

In Figure 3c, we illustrate the magnetic-field evolution 
of four non-degenerate valley-orbit levels by means of 
an elementary model. Each level splits into spin-up and 
spin-down levels in finite magnetic field. We assume 
that the exchange interaction is small in comparison 
to the level spacing. The level crossings that follow 
from our model fit the kinks observed in the first four 
Coulomb peaks. The observed kink positions yield three 
valley-orbit levels which are 0.10, 0.23 and 0.29 meV 
above the lowest ground state level. The extracted level 
spacings for the first four valley-orbit states are then: 
E V0 2 ~ E vol = 0.10 meV; E vos - E V02 = 0.13 meV; 
and £Vo4 — ^vo3 = 0.06 meV. 

5. Spin filling of electrons 5—12 

Finally in Figure 4, we plot the addition energies 
E&dd(N) as a function of B for electrons N = 5 to 
12. Once again, we predominantly observe slopes of 
dE&dd(N)/dB = 0, as expected from Equation 

(2). Occasionally, e.g. at N = 6 7, a segment has 
a slope of ±2g{iB, because the total spin on the dot 
changes by more than ^. This can occur due to many- 
body interactions on the dot and lead to spin block- 
ade [29]. The latter phenomenon could also explain the 
suppression of current in the fifth charge transition at 

b = 2-5 t planum]. 

Also, the picture of alternating spin filling below 
0.8 T no longer holds for N > 4. Unexpectedly, the 
fifth electron is spin up at low magnetic field, while the 
lowest-energy configuration predicts a spin-down state. 
This anomaly could be explained by an extra electron 
in a dot nearby, which alters the spin configuration of 
the main dot. Such a small dot can be created at high 
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FIG. 4: Addition energies of the 5 th to 12 th electron versus 
magnetic field, offset for clarity. Black numbers -2, -1, 0, +1 
and +2 correspond to the slopes of the addition energy in mul- 
tiples of g\±B and reveal the sign of the added spin (indicated 
as red and green arrows) according to equation (2). Coulomb 
peak spacings in gate voltage are converted to energies using 
the lever arm ap, which vary from 0.11 to 0.068 eV/V with in- 
creasing electron number. The measurement was taken along 
the dashed line marked in Fig 3(a) at Vb2 =0.738 V for the 

electron were taken at Vb2 =0.710 V. 



th to 7 th electron. The addition energies of the 8 th to 12 th 



plunger gate voltages, where the potential well differs 
from a perfect parabola. As more electrons are added 
to the main dot, the wavefunctions extend further and 
would have more opportunity to spin-couple to the 
unintentional dot nearby, thus affecting the spin filling 
of electrons. 



6. Conclusion 

The results here show that silicon MOS quantum dots 
can be fabricated with the low levels of disorder neces- 
sary to form well-defined electron spin qubits in a host 
material that can be made almost free of nuclear spins. 
The excellent charge stability allows the spin states of the 
dot to be mapped up to N = 12 electrons and a valley- 
orbit splitting of 0.10 meV to be extracted. A recent 
theoretical study [15] has shown that a valley splitting of 
0.1 meV is sufficient for the operation of a silicon double 
quantum dot as a singlet-triplet qubit, in analogy with 
recent experiments in GaAs [30]. Given that the valley- 
orbit splitting is strongly dependent on the interfacial 
electric field, it should be possible to further increase 
the splitting via appropriate device engineering. Our re- 
sults therefore provide real promise for the realization 
of low-decoherence spin qubits based upon silicon MOS 
technology. 
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